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DISCUSSION OF TESTS RESULTS IN THE DESIGN OF LAMINAR 
AIRFOILS FOR COMPETITION GLIDERS 

Jerzy Ostrowski, Stanislaw Skrzynski and Mieczyslaw Litwinczyk 

Optimization of flight profiles, especially for swift planes, /105 
has usually been achieved from the point of view of reducing lag co- 

'Mi' 

efficient. The carrying power' of such profiles has been augmented 
by appropriate mechanization of the wing in an attempt to attain co- 

I M H . ; i-\ ! c • /• : , . , e.y, l | y» •/. 

efficients of maximal carrying power 1 considerably higher than those 
of the basic profile. 

In glider construction there is a limit on wing mechanization: 
either flap and brake (open class) or brake alone (standard class) . 

This is the reason for dissimilarity of standards set for glider pro- 
files, especially in the standard class. Here the criterion of mini- 
mal lag is still in force, yet the requirements for the size of c 

Z Ttl c9 X 

are on the increase. It is particularly important to achieve high 
c values for small (in the scale of aviation applications) Re 

numbers, starting with values less thin one million. 

For a long time glider builders used NACA laminar profiles. In 
the early 1960's profiles of the Fx (Wortman) type appeared with bet- 
ter properties than the NACA. However there was still a need for com- 
plete documentation of these profiles as well as of a good deal of in- 
formation necessary for constructors' estimates of the features of the 
glider being planned. 

This was the situation when during that decade a joint commission 
came from the Bielsk Glider Factory and the former Aerodynamics De- 
partment of the Warsaw Polytechnic to undertake research aimed at se- 
curing the missing information in respect to the properties of pro- 
files in the range of small Reynolds numbers and to run tests on the 
modifications of such profiles for improving their features. For this 


* 

Numbers in the margin indicate pagination in the foreign text. 

2P 


1 




Fig. 1 


purpose we built a small la- 
minar tunnel, providing for 
its regulation, and we planned 
and set up the instrumentation 
(Fig. 1) . We then carried out 
preliminary flow characteris- 
tic tests, worked out a meth- 
od of measurements, deter- /106 
mined improvements on the 
tunnel and assessed a number 
of factors affecting the range 
and results of measurements. 



Fig. 2 

Key: a. turbulence 

b. separation 

c. adherence 



Fig. 3 


Key: a. theoretical 
b. actual 


Analysis of the reflux for various pro- 
files and the existing experimental data (e- 
ven though at the beginning of the testing 
these were very scant for small Re numbers) 
indicated the fairly real influence of cer- 
tain flow phenomena in the laminar border 
layer, particularly in the zone of transition 
between laminar and turbulent flow, the val- 
ue for c , the course of separation and 

z max 

the plotting of the c and c curves. These 

z x 

flow phenomena, until now studied without 
much system, occurred in the shape of fully 
developed "laminar bubbles". 

As a function of profile formation. 

Re number and surface unevenness these 
bubbles arise in the area of narrower or wi- 
der angles of friction and occasion either 

changes in the course of the c and c curves 

z x 

or else a slow or rapid separation process. 
Fig. 2 shows how a laminar bubble occurs. 

It is characterized by separation of the la- 
minar border layer and adherence of this 
layer as a turbulent layer. In its fully 
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Fig. 6 


developed form the bubble be- 
comes a separation area with 
considerable recirculatory 
flow. 

The term "bubble" does 
not do justice to the charac- 
ter of the observed phenome- 
na. Although the laminar bub- 
bles form in the area directly 
preceding separation in all /107 
test cases, we observed a num- 
ber of deformation phenomena 
in the velocity field without 
recirculatory flow. In a real 
way this deformation altered 
flow structure in the layer, 
which may be considered as the 
earliest stage of bubble evolu- 
tion. Then too the size and 
shape of the bubbles themselves 
were different. In measuring 
test bubbles on profiles from 
an 0.3 m chord to aim chord 
we were able to differentiate 
bubbles several hundredths of 
a millimeter thick and on the 
order of 1 mm. Their length 
varied from several to several 
hundred mm. The term long bub- 
bles is used if the length is 
from several dozen to several 
hundred times their thickness 
[2] . Short bubbles usually 
occur near the nose of the pro- 
file and represent local sepa- 
ration in miniature, such as 
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for example one notices when there is 
separation in a sharp edge environ- 
ment. Such separations of a local na- 
ture are included in this classifi- 
cation because their thickness is much 
greater than that of a border layer 
(often man times thicker) . 


Fig. 7 



Fig. 9 


A particularly large amount of 
time was devoted to the testing of la- 
minar bobbles , since the elaboration of 
appropriate forms for laminar airfoils 
results in the selection of a form 
where further lamination retards the 
appearance of bubbles and reduces them 
in size. 

Experimental elaboration was un- 
dertaken from the viewpoint of the 
fact that, despite appearances, there 
existed a great deal of imperfection 
in profile theory, which, although it 
has become the formally correct scheme 
and the one enshrined in descriptions 
of flow phenomena according to the 
classical concept, yields such diver- 
gent results in respect to values got- 
ten experimentally, that neither the 
pressure distribution derived from a 
theoretical description nor the distri- 
bution of tangential stresses may be 
used as the basis for an assessment 
of the effect of changes in airfoil 
contours on their aerodynamic charac- 
teristics. Here there is especially 
a great divergence in the Slope of 

the actual curve c and of the theor- 

z 
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Fig. 11 


MCA 65-212 



Fig. 12 

etical curve (Fig. 3) , particularly great for uneven profiles for 
which the value dc /da obtained theoretically increases greatly with 

Z 

profile uneveness. In actuality this increase is much less and thus 

there is greater divergence. Further we note a difference in the 

case of asymmetric profiles between zero angle values when carrying 

power is a _ Q . For the profile tested in real fluid the zero angle 
z 
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Fig. 14 


of the carrying power is less than one would theoretically expect. 

Tests for a way out of the impasse by the introduction, as a /108 

reference value, of the coefficient c in place of the friction 

z 

angle a [5] and transformation of theoretical calculations to suit 
real conditions are not justified in view of the divergence found 
in the results. In the first place, the angle of friction was at 
the time undefined and in the second place, with such a great dif- 
ference A a the theoretical velocity field will be quite different 
from the actual one and cannot be used as the basis for an analysis 
of flow phenomena in a border layer. In the test area for these phe- 
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nomena the theory of the lami- 
nar layer then no longer yields 
any information in view of the 
assumption of pressure stabili- 
ty in a normal direction and 
the gaps in the turbulence the- 
ory, particularly in the area 
of processes that occur in the 
regions where the stability of 
the laminar layer is lost. 

Fig. 4 presents diagramma- 
tically a typical curve for la- 
minar lag.. Here one may dis- 
tinguish five distinct sections 
characteristic of profile flows. 
Section 1 corresponds to reflux 
where a large portion of the up- 
per surface experiences laminar 
flow in a layer whereas the en- 
tire lower surface is marked by 
turbulent flow; section 3 is 
where laminar flow occurs over 
a large part of the upper and 
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lower portions of the profile, while section 5 is the opposite of 
section 1. In sections 2 and 4 we find more or less rapid changes 
in the value for the coefficient c due to corresponding lamination 
of flow on the lower surface or flow disturbance on the upper. The 

size of such sections and the course of the a curve in them is a 

x 

function of the rapidity with which the transition points shift (or 
more exactly transition zones) with changes in the angle of friction. 

The shifting process for transition zones is a function of a number 
of factors, precisely of the kind and formation site of laminar bub- 
bles as well as of velocity field deformation in the border layer. 

They appear on the lower surface with small and negative c values 

z 

and on the upper portion of the profile when those values are large. 

For an understanding of the reasons for the appearance of such 
phenomena it seems useful to discuss the process of laminar bubble 
formation. Here we may make use of the scheme presented in [3] il- /109 
lustrating the occurrence of a bubble on a circular profile in the 
area of critical Re numbers. Fig. 5 shows how this scheme is ap- 
plied to a flight profile. When the Re numbers are small and sub- 
critical, laminar separation occurs near the site of a change in the 
symbol of the pressure gradient and this change is a function of the 
shape and arrangement of the profile. The separated stream is sub- 
jected to turbulence at a certain distance behind the airfoil. As 
the Re number increases the turbulent region of the separated stream 
approaches the surface of the airfoil (Fig. 5b) , whereas if it ap- 
proaches the surface at a correspondingly short distance the flow 
is shut off (Fig. 5c) as the result of turbulent diffusion and cre- 
ates the characteristic laminar bubble shown in Fig. 2. With further 
increase in the Re number, depending on the shape of the airfoil con- 
tour (change in pressure gradient) the bubble may disappear with the 
appearance of a transition zone from laminar to turbulent flow in 
the border layer without separation. In some cases one may note 
considerable deformation of the velocity field in the transition zone. 

In other cases, when pressure gradients are quite small, we observe 
no velocity field deformation in the layer. There may then exist, 
depending upon Re number and contour formation, different size bub- 
bles and various more or less obvious deformations in the border lay- 
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er. From the change in the friction angle of the profile the ob- 
served deformations may pass down the contour assuming different 
forms. In some cases large bubbles may appear at fairly narrow an- 
gles of friction and persist over a very wide range of Re numbers; 
in other cases we may find quite small bubbles among the borders of 
friction only with large angles of friction. Figures 6, 7 and 8 pre- 
sent bubbles in order: a bubble at the margin of friction, a bubble 
near the point of maximal thickness (in the diffusion zone) as well 
as a broad area comprising deformation of the laminar layer. Bubbles 
are likewise observed on wings with full spread and different con- 
tours . Figures 9 and 10 show in succession a short bubble at the 
margin of friction on a wing tilted to y = 30° and a long bubble on 
one tilted to y = 60°. 

Their effect however on separation phenomena is different, less 
in respect to the tridimensional nature of reflux from a tilted wing. 

/HI 

Fig. 11 gives a schematic representation for two classical cases: 
a thick profile with a fully developed bubble behind the point of ma- 
ximal thickness and a slender profile with a small bubble at the nose. 

On certain types of profile we may observe, even with narrow an- 
gles of friction, flow deformation in the laminar layer comprising re- 
moval of the layer from the wall of the profile and the setting up of 
a very thin sublayer with a small linear velocity gradient. Such de- 
formation does not bring about the loss of layer stability. The lay- 
er behind the deformation area still retains a laminar velocity pro- 
file. Such a case of deformation may be regarded as the very earli- 
est stage of a laminar bubble, since it may be seen at those sites on 
the contour where, when the friction angle increases, we observe a la- 
minar bubble or deformation that occasions flow turbulence. 

It follows from the above description that over a wide range of 
changes in the geometry of the profile and Re numbers we may note a 
variety of laminar bubble forms, both quantitatively and qualitative- 
ly, and deformation of the velocity field in the layer. They are il- 
lustrated in the Figures presented. Figure 12 shows successively 
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the forms of evolutLonof a bubble with a change in the angle of fric- 
tion and illustrated by means of lines of steady velocity. Fig. 12 
(left) shows field deformation without local separation (profile 
NACA 65-212, friction angle a = 6) . Fig. 12 (middle) presents in 
turn an already fully developed laminar bubble (a = 7.5°) and Fig. 

12 (right) the same bubble for a friction angle a = 9°. 

Fig. 13 illustrates velocity field deformation in the setting of 
a laminar bubble and likewise involving the transition zone from la- 
minar to turbulent motion; Fig. 14 shows the transition zone in the 
deformation area (without separation) . 

Fig. 15 represents a velocity field in the transition zone with /112 
no separation resulting from the great changes in pressure gradient. 

Bubbles may likewise be visualized by covering the area surface 
with a matching adhesive fluid. No matter how much the presence of 
the liquid changes local reflux conditions, the sites of a bubble com- 
ing into being may be determined and also its size. However visuali- 
zation does not permit us to determine whether we are dealing with a 
developed bubble or with field deformation. In the latter case a con- 
siderable reduction in tangential tensions favors accumulation of 
oil and this may result in border layer separation at the site. 

As we see from the phenomena discussed previously, laminar bub- 
ble formation constitutes a phenomenon directly preceding separation. 
Hence the basic question in planning a good profile has been the se- /113 
lection of profile shape so as to retard the appearance of such bub- 
bles or to preclude their development. 

A good glider profile should have these features: the least pos- 

sible lag in areas of high velocity (as most of the contour is bathed 

in a laminar layer) ; the greatest possible c and in addition 

z max 

the beginning of the saddle shifted far toward the side of negative 
angles of friction, so that at very high velocities (small coeffi- 
cients of carrying power, such as c ^ 0.1) the lag coefficient c is 

z x 

located on the borders of the saddle. Formation of a long laminar 


10 



1 


layer requires shifting the maximum width of the profile far to the 
rear and this in turn suggests the necessity of providing a contour 
of great curvature in the diffusion zone, so much the more since in 
consideration of the need for obtaining a high c ^ value the 

shaft of the profile's skeletal line must be considerable (3-4%) . 

Sharp curves in the contour in the diffusion zone favor the forma- 
tion of bubbles and these cause partial separation even at relative- 
ly narrow angles of friction (small c max ) • 

Shifting the position of the saddle toward the smallest possible 
c' z values makes profile asymmetry difficult (a^ “ -4) and this, in 

the face of the position of the accumulation point above the border 
of friction in this area of friction angles, creates great diffi- 
culty in constraining reflux from the nose without separation or de- 
formation that cause turbulence of the layer on the bottom surface 
of the contour. 

We see from this that all the changes of contour that have a fa- 
vorable effect on the c . value, the coefficient for the c moment 

x mm m 

and the extreme position of the saddle have an adverse effect on the 

size of c and the course of separation. 

z max 

Reconciliation of these discrepancies in an optimal way would re- 
quire study of profiles of different shapes and assessment of the ef- 
fect of geometrical differences for individual portions of the con- 
tour on pressure distribution, flow phenomena and actually on its 
characteristics . 

We used here the method of successive changes in the shape of se- 
lect portions of the contour while retaining unaltered basic geomet- 
ric parameters of the test profiles. In this way, through a series /114 
of modifications, we achieved a profile with definite aerodynamic 
characteristics. Then, altering step by step the basic geometric 
parameters (e. g. position of maximum thickness) and the contour, we 
obtained profiles whose characteristics, along with the results of 
reflux testing, put us in a position to assess the effect of changes 
in geometric parameters and the shapes of test profiles on their aer- 
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odynamic characteristics. Profiles obtained in the course of such 
tests and exhibiting certain characteristics in relation to previ- 
ously tested ones were described and numbered so as to produce a se- 
ries of profiles NN comprising 10 to 20 and including several result- 
ing from modifications of NACA profiles and those of the F^ series. 

As a result of this testing program we obtained information in 
regard to the possibility of setting up pressure distribution in a 
manner that would provide stability in the marginal layer over a wide 
range of friction angles (retardation of occurrence and diminution 
of occurring deformation of laminar flow) . 

It appeared that a profile meeting the above requirements had to 
incorporate -- in contrast to classical hypotheses (contour selec- 
tion for a required pressure distribution at a given friction angle) 
— definite changes in pressure distribution along a chord with a 
change in the angle of friction. 


Without entering into the specifics of individual stages we will 
discuss below three types of characteristics with corresponding chan- 
ges in pressure distribution and forms of deformation of the layer 
that have a bearing on individual characteristic types. 

Fig. 16 presents characteristics of the first type. As the Re 
number goes down separation shifts toward wider friction angles 
(curves 1, 2 and 3) and even reaches the area of the laminar saddle. 
There is no displacement here of the transition zone in the border 
layer toward the front with mounting Re numbers. Deformation of the 
velocity field, which in this case occurs in the rear portion of the 
profile behind the point of maximal thickness, changes to a long lami- 
nar bubble as the angle of friction widens. The c characteristics /115 

z 

for this type of profile present two maxima in the range of small Re 
numbers: the first, when there is separation in the diffusion zone, 

caused by the bubble discussed above; the second occurring with 
wide friction angles following the creation of a second (short) bub- 
ble at the profile nose. Typical representatives of this type of 
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profile are the thick profiles NACA of the 65 and 66 series. 

If there is nothing altered about the profile shape, so that in 
the rear portion only a weak bubble is formed, then as the Re number 
changes there will be a change in the c characteristic in the man- 
ner shown in the Figure. A second maximum for this characteristic 
appears and at the corresponding Re number the bubble in the diffu- 
sion zone disappears. We then get the c z characteristic shown in 
Fig. 16 by the heavy line (shown likewise in Fig. 4) . The second 
maximum for this characteristic enters the area of usable angles of 
friction. In the area where the bubble appeared we note weak deform- 
ation of the laminar velocity field and, for angles corresponding to 
to inception angle of the bubble at lower Re numbers, a break or leap 
in the c angle. When the Re numbers are very small, ^ 0.5 *10 6 , the 

Z 

form of the lag curve deviates from what has been discussed — we 

note a second c minimum with wide angles of friction. In the region 

x 

of this minimum there is a large bubble (curve 5) . 

Fig. 17 shows changes in pressure distribution with an increase 
in the angle of friction for this type of profile. We see clearly 
that as this angle widens there is a rapid increase of hypopressure 
around the nose that creates conditions favoring the appearance of a 
very long bubble which produces instant separation in the rear part 
of the profile. This separation results in a reduction of carrying 
power and therefore also of circulation. There is a change in the 
position of the point of entry and in wing reflux. In this context 
flow is sucked against the forepart of the profile. Hypopressure 
in this area increases with a further increase in the angle of fric- 
tion, creating a distribution typical of a bubble arising near the 
margin of friction (typical flattening out behind the point of max- 
imal vacuum) . This is shown on the diagram with heavy lines. 

Fig. 18 shows the characteristics of profiles where there is no 
formation of a long bubble of the type discussed. Here we note a /116 
gradual increase of lag from the end of the saddle and separation 
occurs as the result of the bubble formation at the friction border. 
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In the area of friction angles, where 
we observe increasing lag the transi- 
tion zone on the upper portion of the 
profile shifts forward. Here changes 
in pressure distribution as the angle 
of friction widens are characterized 
by a much evener increase of vacuum 
and a forward shift of the hypopressure 
maximum of Fig. 19. The bubble forms 
at considerably greater angles of fric- 
tion so that c is much greater. 

This is the type of characteristic we 
find in profiles such as F^61-163, 
F^63-167, NN 6 and thick asymmetric 
NACA profiles of the 64 series. 

A third type of characteristics 
(Fig. 20) , the least studied because 
of the difficulty of measuring very 
small bubbles that arise near the nose, 
differs from the preceding in the gra- 
dual shift of the hypopressure maximum 
towards the friction border (compare 
Figures 17 and 19 with 21) . Admitted- 
ly this leads to a rapid monotone in- 
crease in lag almost from the c 

x min 

(here there is no typical saddle form- 
ation) but formed bubbles do not ap- 
pear in this case, since changes in pressure gradients are less over 
a wide range of friction angles than in the previous cases. Here a 
special feature is likewise constraint in the appearance of a small 
laminar bubble at the border of friction even in the area of negative 
pressure gradients. Such a bubble disturbs the layer, making it sta- 
tic and retarding its separation and the latter occurs only when, with 
an increase in the angles of friction, the bubble shifts toward the 
side of positive pressure gradients. This phenomenon is accompanied 
by changes in pressure distribution in the environment of the borders 
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of the nose as shown in Fig. 21. Both 
these factors — constraint of the cha- 
racteristic changes in pressure distri- 
bution with a change in the angle of 
friction and constraint of the very 
small bubble referred to in the area 
of negative pressure gradients — here 
are decisive in respect to the large 
value of c achieved for this type 

of profile. The real question in the 
attainment of correspondingly low lag 
is the constraint of reflux fulfils /117 
the requirements of the above condi- 
tions when the maximal thickness of 
the profile is shifted as far back as 
possible. 



Profiles that are still laminar 
and possess this sort of characteristic 

are NN7, NN8 and Wortmann F 63-168 with 

x 

a slight modification of the nose in 
accord with test results (profile NN8 
is a modification of F^63-168 compris- 
ing the forepart and the flow portion, 
due to which we achieved moment reduc- 
tion by about 30%, better placement of 

the saddle and a c . less than over 

x mm 

5%) . A comparison of characteristics 
for the last two profiles is given in 
Fig. 22. 


Fig. 21 


For a general assessment of the 
aerodynamic properties of profiles re- 
presenting types of characteristics 
the following data may be helpful: first type — profiles optimized 
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Fig. 24 

Key: a. testing by the Stuttgart Institute 

testing in Jasien 


from the viewpoint of the minimum value, profiles which present at 

Re 'v 1.5xl0 6 : c . ^0.0055 and c 'v 0.9; the second type for 

x mm z max 

the same Re number presents: c . ^ 0.006 and c -v 1.3; the 

x mm z max 


third type c . ^0.0065 and c ^1.45, 

x mm z max 


A comparison shows that when there is a relatively small, about 

15%, increase in the c . n . ...... 

x mm an approximately 60% increase is obtained 

for c 

z max . 

Evaluating the results on the basis of practical use we must pay 
attention to the very high sensitivity of the test profiles in respect 
to changes of shape in areas where large contour curves appear and in 
particular the nose and its environs. Very slight inaccuracies in pro- 
jection of a contour and even an insignificant amount of surface rough- 
ness may occasion considerable changes in the formation of the lami- /119 

nar saddle as well as an adverse increase in c (up to 10%) and a 

x ihin 

lowering of c z max . Inaccuracies in producing the model are the chief 
reason for observed differences in measurement results on the model of 
this very profile. Fig. 23 shows by way of example charts drawn up 
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for two models of the profile of the 13% NN 8-13. Fig. 24 presents 
results of testing done on airfoil NN 7MK in the wind tunnel of the 
Stuttgart Institute on which we have superimposed measurement re- 
sults of our own for this profile. Here we see clear deterioration 
of the laminar saddle (it appears in the curves plotted for our own 
measurements) , although on the other hand when we repeated our own 

measurements we also got the c curve without formation of the sad- 

x 

die. One can also note good agreement with measurement results from 
countries outside Poland, e. g. for the Wortmann F^61-168, in respect 
to those obtained in our own wind tunnel. Airfoil sensitivity suggests 
the need for great precision in determining the geometry of the mod- 
el and transfering it to a real wing while maintaining very great ac- 
curacy. 
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